Mass Models for Spiral Galaxies from 2-D Velocity Maps 
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ABSTRACT 

We model the mass distributions of 40 high surface brightness spiral galaxies inside their 
optical radii, deriving parameters of mass models by matching the predicted velocities to observed 
velocity maps. We use constant mass-to-light disk and bulge models, and we have tried fits with 
no halo and with three different halo density profiles. The data require a halo in most, but not 
all, cases, while in others the best fit occurs with negligible mass in the luminous component, 
which we regard as unphysical. All three adopted halo profiles lead to fits of about the same 
quality, and our data therefore do not constrain the functional form of the halo profile. The 
halo parameters display large degeneracies for two of the three adopted halo functions, but the 
separate luminous and dark masses are better constrained. However, the fitted disk and halo 
masses vary substantially between the adopted halo models, indicating that even high quality 
2-D optical velocity maps do not provide significant constraints on the dark matter content of 
a galaxy. We demonstrate that data from longslit observations are likely to provide still weaker 
constraints. We conclude that additional information is needed in order to constrain the separate 
disk and halo masses in a galaxy. 



Subject headings: galaxies: fundamental parameters - 
matter — galaxies:stellar content 

1. Introduction 

The mass of a galactic stellar population has 
proved to be a particularly difficult quantity to 
determine. In principle, an accurate value would 
place a constraint on the star formation history, 
stellar mass function, and metallicity. But it is 
perhaps more important for dynamics because it 
affects estimates of the dark matter content of the 
inner parts of galaxies. There are two possible ap- 
proaches to estimating the stellar mass: dynami- 
cal measurements and comparison of the observed 
broad-band colors or spectra with theoretical pre- 
dictions from population synthesis. 

Theoretical models to synthesize a stellar pop- 
ulation {e.g., Arimoto & Yoshii 1986; Bruzual & 
Chariot 1993; Worthey 1994; BeU et at 2003) use 
established stellar models and assumptions about 



- galaxies:kinematics and dynamics — galaxies:dark 

the likely mass function, star formation history, 
etc. to predict the luminosity in a broad-band 
color. The expected luminosity in almost any color 
band is insensitive to the low-mass slope of the 
adopted mass function, since large proportions of 
low-mass stars can be added or removed with lit- 
tle effect on the total luminosity. The resulting 
mass-to-light values, T, are thought to be reliable 
to within a factor of 2 (BeU & de Jong 2001). 

The dynamical approach is scarcely more suc- 
cessful. The typical circular flow pattern of gas in 
a spiral galaxy affords a rather precise estimate of 
the overall mass profile within the radial range ex- 
plored by the data. But it is much harder to deter- 
mine the fraction of the enclosed mass that comes 
from the starlight, since an unknown fraction of 
dark matter is also expected to contribute. A com- 
mon simplifying assumption is that the disk mass 
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profile traces that of the disk light with a constant, 
but unknown, T. But van Albada et al. (1985) and 
Lake & Feinswog (1989) stress that, for an ax- 
isymmetric flow pattern, the disk and halo masses 
are almost completely degenerate, even for galax- 
ies with the best available rotation curve data from 
neutral hydrogen observations. 

An additional factor that contributes to this dif- 
ficulty is the unknown density profile of the dark 
matter halo. Previous studies have tried a num- 
ber of functional forms, some motivated by cos- 
mological simulations, but with no decisive con- 
clusions {e.g., Navarro 1998; .Jimenez, Verde & Oh 
2002). It is scarcely surprising that very weak con- 
straints result from fitting a 1-D rotation curve as 
the sum of a contribution from the luminous mat- 
ter, with an unknown mass scaling, plus a contri- 
bution from dark matter with both an unknown 
functional form and an unknown density normal- 
ization. 

Weiner, Sellwood & WiUiams (2001b), Weiner 
(2003), and Kranz, Slyz & Rix (2003) were able 
to obtain more precise estimates (within ^ 10%, 
Weiner et al. 2001b) of the stellar disk mass 
through modeling the non-axisymmetric fiow pat- 
tern in disk galaxies. These were time-consuming 
studies, requiring a large grid of hydrodynamical 
simulations over the two unknowns, T and the pat- 
tern speed of the bar or spiral. Alternatively, Bot- 
tema (1997) measures the velocity dispersion of 
the stars in the disk, to obtain a disk mass esti- 
mate (see also Vcirhciijen et al. 2003). While these 
more challenging methods may yet be needed to 
break the disk-halo mass degeneracy, we here ex- 
amine the question of whether fitting simple ax- 
isymmetric models to high-q\iality, 2-D optical ve- 
locity maps offers any tighter constraints over pre- 
vious attempts at fitting 1-D rotation curves from 
longslit data or lower-resolution 2-D maps from 
aperture synthesis. 

Velocity maps at low spatial resolution have 
long been available from neutral hydrogen aper- 
ture synthesis observations {e.g., Bosma, van der 
Hulst & Sullivan 1977). More recently, velocity 
maps have become available at very high spatial 
resolution from mm-wave observations of molecu- 
lar gas (Jogee et al. 2002; Sofue et al. 2003; Si- 
mon et al. 2003). Most optical data has contin- 
ued to be obtained using longslit spectra {e.g., 
Courteau 1997; Sofue & Rubin 2001; Vogt et al. 



2004), but a few maps are available from integral 
field units; Corradi et al. (1991) and Palunas & 
Williams (2000) used Fabry- Perot maps, Andersen 
& Bershady (2002) have used fiber bundles, and 
the lens array spectrograph SAURON has been 
used to investigate early type galaxies (Davies et 
al. 2001). 

Palunas & Williams present the largest sample 
of galaxies with homogeneously gathered Fabry- 
Perot maps and /-band photometry available so 
far. They showed that the inner rotation curves 
could be fitted by constant Tj mass models and no 
dark matter, but with a rather wide range of T/ 
values, when mass discrepancies in the outer parts 
were excluded from the fit. We here present a re- 
analysis of these galaxies to determine what can be 
concluded about the likely dark matter content of 
galaxies from fitting axisymmetric models to these 
data. 

Our objective here is to try to determine 
whether our data place meaningful constraints 
on the dark matter content of the galaxies in our 
sample and whether one halo profile or another 
yields a sigificantly better fit within the radial 
range spanned by our data. These objectives are 
quite different from testing theoretical predictions 
from some model of galaxy formation. 

We find that 2-D velocity maps do provide some 
advantage over 1-D rotation curves, in that they 
reduce the statistical uncertainties in estimated 
parameters. However, the disk mass, and there- 
fore also the dark matter content of a galaxy, is 
still strongly dependent on the adopted halo den- 
sity profile. We support these conclusions by fit- 
ting mass models to both the 2-D velocity map 
(§5) and to a 1-D rotation curve (§7) which we 
derive from data along only the major axis of our 
velocity maps. 

2. Data 

The data used in this paper were kindly pro- 
vided by Povilas Palunas and Ted Williams (see 
Palunas & Williams 2000, hereafter PW, for 
a complete description of the data and reduc- 
tions). We have /-band images and Fabry-PcTot 
(hereafter FP) velocity maps for 74 high surface- 
brightness spiral galaxies. PW report the aver- 
age seeing to be « 1.5". We adopt the mag- 
nitudes, luminosities, and distances reported by 
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PW; distances are simply Hubble distances as- 
suming Ho = 75 km s~^ Mpc~^. 

Galaxies in this sample were selected to lie in 
a particular part of the sky, to be of later Hubble 
type, to not be classifed as barred in the RC3 cat- 
alog (de Vaucouleurs et al. 1991) (although many 
are, in fact, barred), to not be obviously disturbed, 
to have a suitable angular size for the FP instru- 
ment, and to have redshifts in the range 4000 - 
5000 km/s. The data were taken to obtain optical 
Tully-Fisher distances of galaxies to map depar- 
tures from the Hubble flow in the general direction 
of the Great Attractor, which may have resulted in 
a somewhat larger than normal fraction of galaxies 
in clusters, although roughly half (35 of 74) are in 
the field. See § 4.1 for a discussion of the impact 
of the Great Attractor on the distances adopted 
for these galaxies. Apart from this, the sample 
seems unlikely to have any sinister biases relevant 
to their mass distributions. 

We have created mass models for a subsam- 
ple of 40 galaxies. The T- types of this subsam- 
ple (from the RC3) range between 2 (Sab) and 10 
(Irr) with most being 4-6 (Sbc-Sc). The radii at 
which the surface brightness is 23.5 mag/arcsec^ 
-R23.5 range from 24 to 104". Their average lumi- 
nosity is 3 X 10^'^ Lq and ranges between 8 x 10^ 
and 8 x 10^° Lq. The average I- band central sur- 
face brightness is 19.6 /i, with the brightest galaxy 
at 18.3 ^ and the faintest at 21.3 /it. Maximum 
circular speeds for these galaxies range from ap- 
proximately 100 to 300 km/s with an average of 
170 km/s. The 34 galaxies in their sample which 
we have omitted from this study are listed in Ta- 
ble 1. The majority were discarded because they 
are close to edge-on, having inclinations greater 
than 75°, which raises a number of difficulties: (1) 
The finite thickness of their disks becomes a sig- 
nificant fraction of the projected minor-axis and 
invalidates our simple razor-thin disk model (see 
§3.1). (2) Extinction problems become severe, and 
dust corrections large and uncertain. (3) The ve- 
locity map is less well resolved, and gas emission 
on the minor axis can more easily skew the esti- 
mated velocity from the major axis, which could 
lead to an estimated rotation curve that may be 
more slowly rising than is really the case. A few 
galaxies were discarded for other reasons, such as 
a lack of kinematic data in their central regions. 

We did not discard galaxies on the basis of their 



morphology and have persisted in fitting axisym- 
metric models to a few that manifestly have strong 
bars. The reason for this is that almost all galax- 
ies in the sample have weak non-axisymmetric fea- 
tures, and judgements of which galaxies to discard 
would have been arbitrary and subjective. Model- 
ing of the non-axisymmetric flow pattern would be 
superior, but it is a major task for even one galaxy 
(Weiner, Sellwood & Williams 2001b; Kranz, Slyz 
& Rix 2003) which we do not attempt here. 

Throughout this paper, we fix the position an- 
gle 4>0) inclination angle i, and systemic velocity 
Vs of each galaxy at the values deduced from ear- 
lier fits of axisymmetric flow patterns to our kine- 
matic data (Barnes & Sellwood 2003). Briefly, we 
found the values of these parameters that mini- 
mized between the 2-D velocity map from the 
FP data and a planar axisymmetric flow pattern. 
The fit therefore also yields a set of mean orbital 
speeds, and their uncertainties, at equally-spaced 
radii which we show as the "data" in subsequent 
figures. Since we fit axisymmetric mass models, 
we could, in principle, fit this 1-D rotation curve, 
which is basically the speed averaged in annular 
bins, provided we compute the weights correctly. 
However, our use of the biweight (§3.3) requires a 
fit to the full 2-D velocity map. 

The use of kinematically-derived projection pa- 
rameters is possible only because we have velocity 
maps. Rotation curves derived from long-slit ob- 
servations must adopt inclination and major-axis 
position angles from the brightness distribution in 
the galaxy. We discuss how much this might affect 
the results in §7. 

3. Methods 

We build a mass model for each galaxy in our 
sample based on the measured light distribution in 
the /-band and a possible dark matter component, 
and then compare the 2-D velocity map expected 
from gas on circular orbits in our model with that 
observed from the Ha emission. 

We assume separate constant /-band mass-to- 
light, T/, for the disk and bulge components and 
neglect the possibility of T/ gradients within each 
component. This is the usual assumption, despite 
established color gradients in galaxy disks {e.g., de 
Jong 1996). We justify this assumption because 
the unknown correction for T/ gradients will have 
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a small effect, as shown by Weiner et al. (2001a): 
the rotation curve shape for the disk is largely de- 
termined by the mass distribution in the dense 
inner parts, and changes to the density gradient 
at larger radii have a minor effect. 

We also neglect any contribution from gas. In 
previous models where HI maps are available {e.g., 
Broeils 1992; Verheijen 1997), the contribution 
from atomic gas to the rotation curve is generally 
very small in the inner parts, but rises toward the 
optical edge - to < 10% of the central attraction in 
the most extreme case. Molecular gas is generally 
more concentrated towards the inner parts {e.g., 
Regan et al 2001) where its contribution may be 
somewhat larger, but we could not find any infor- 
mation on the molecular gas content of our galax- 
ies. HI fluxes for 28 galaxies in our sample are 
available from NED.^ Since we find gas mass frac- 
tions typical of the galaxies in these other stud- 
ies, neglect of the gas contribution should result 
in only a slight systematic overestimate of the T/ 
for the stellar component, unless the dark matter 
fraction is large. 

We consider three different radial density pro- 
files for the dark matter halo, which we assume to 
be spherical for simplicity. Two of our halo func- 
tions are standard: the pseudo-isothermal model 
with a constant-density core and a model with a 
broken power law; the third is a simple, scale- 
free power-law. Our objective is not to test any 
particular model of dark matter halo formation, 
but to determine how well our data constrain the 
halo contribution to the rotation curve in the disk 
plane. It is possible that halos are spheroidal, or 
even triaxial; the unknown shape of the halo may 
alter the inferred moan density profile slightly, but 
our data cannot constrain such subtleties. The 
calculation of the contribution to the circular ve- 
locity from a spherical halo model is straightfor- 
ward. 

3.1. Bulge-Disk Decomposition 

Since we wish to allow the disk and the bulge 
to have separate T/ values, we must separate 
the light distribution into these two components. 

iThc NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronau- 
tics and Space Administration. 



While we do this with care, the central attrac- 
tion is not strongly sensitive to the precise decom- 
position, because we ensure that the T/ of the 
two components in our models do not differ by a 
large factor. (We regard it as unlikely on physical 
grounds that Tbj ^ '^d,i', we discard all such 
solutions that emerge from an unconstrained min- 
imization, often with T^,/ 0.1, and repeat the 
fit with the constraint that '^b,i = '^d,i-) We 
also ensure that all the light is included in one 
component or the other. Therefore the total lu- 
minous mass profile remains nearly the same no 
matter how the light is divided. The central at- 
traction depends slightly on whether the mass is 
in a flat disk or a spheroid, with a maximum possi- 
ble difference in circular speeds of ~ 15% (Binncy 
& Tremaine 1987, §2.6), and the fraction of bulge 
light is generally < 10%. Thus the rotation curve 
shape from the luminous matter changes little as 
the division of light between the disk and bulge is 
adjusted. 

Our decomposition method is described in 
Barnes & Sellwood (2003). Briefly, we assume 
the disk to be inflnitesimally thin and axisymmet- 
ric, while the bulge is spheroidal, implying that 
our model isophotes for both components should 
be elliptical but have differing ellipticities, e. We 
assume! the equatorial plane of the bulge is pro- 
jected at the same inclination i and position angle 
00 as the disk plane, which are the values derived 
from modeling the velocity map (§2). We allow 
the bulge to be either oblate or prolate. This 
choice varies from galaxy to galaxy, but most of 
our models are oblate. While we do not assume 
any functional form for the surface brightness pro- 
flle of the disk - it is represented as a table of 
values - we follow PW and model the bulge as 
a sum (up to 3) of Gaussian light distributions. 
Detailed photometry of bulges with high spatial 
resolution {e.g., Balcclls et al. 2003) suggests that 
Gaussians may not be the most appropriate, but 
we find that a sum of three Gaussians can repro- 
duce, e.g.. Scrsic profiles (Sersic 1968) reasonably 
well. Though unorthodox, Gaussians have two ad- 
vantages: the surface brightness is guaranteed to 
decrease faster than any likely disk profile at large 
radii, and they afford a convenient calculation of 
the central attraction (see §3.2). Our assumptions 
arc less restrictive than in schemes that utilize 
specific functional forms for the luminous compo- 
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nents. 

We first attempt to fit the photometric image 
with a disk alone. We then repeat the fit with 
the addition of a single Gaussian bulge component 
and compare the resulting value with its value 
from the disk-only fit. (The dominant source of 
uncertainty is simple photon counting.) If the dif- 
ference is significant, as indicated by an F-test, we 
continue to add a second and possibly a third com- 
ponent, accepting each only if the reduction in ^ 
is large enough that there is a better than 95% 
chance that the extra parameters are required. 
Otu' fitted photometric model typically has final 
reduced (see § 3.3) in the range 1 ^; Xr 

Our best-fit model has the following parame- 
ters; a disk ellipticity eu that is related to the incli- 
nation i of the galaxy through cos i = 1 — e/j, a set 
of disk intensities {/n} at equally-spaced radii, the 
bulge ellipticity es, and the Gaussian central in- 
tensities {/s} and scale lengths {tb) of the bulge 
component (s). 

3.2. Stellar Velocity Field Calculation 

As mentioned earlier, a Gaussian form for the 
bulge simplifies the calculation of the central at- 
traction. Since an axisymmetric Gaussian volume 
light distribution projects into a Gaussian surface 
brightness distribution (Stark 1977), the bulge lu- 
minosity density can be derived from the sum of 
the bulge components through 



j(TO^) 
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where m = i?^ -t- (2/(1 - cb))^) and / = cos^ i -\- 
siT?i/q^g. The apparent bulge axis ratio, qm.Bj 
is related to the bulge ellipticity, es = 1 — qm,B 
and the intrinsic axis ratio of the bulge qb is given 
by ^B — ilm B " '^'^^'^ «)/ sin^ i. Finally, the bulge 
eccentricity es is given by. 
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oblate 
prolate. 



The circular velocity in the equatorial plane due to 
the bulge is therefore (Binney & Tremaine 1987), 



vIb{R) = 47rG^l-e|TB,i^' 



j{rn?)m?dm 



where Tb,i is the mass-to-light ratio. 

We evaluate disk circular velocities using the 
method based on Hankel transforms for a thin disk 
(Binney & Tremaine 1987), 
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where E is the surface mass density, T^j is the 
mass-to-light ratio, and Jq and Ji are cylindrical 
Bessel functions. The cosi factor is an inclination 
correction for the surface brightness, that assumes 
a thin disk with no radial dependence to the ex- 
tinction law. We solve Eq. (3) numerically, using 
a kernel approach. Since we neglect an (unknown) 
thickness correction, which would slightly reduce 
the circular speed for the same surface density, our 
Tjj,! values will be systematically underestimated. 

3.3. Velocity Field Fitting 

We fit our mass model to the galaxy by compar- 
ing the projected velocity map of our axisymmet- 
ric model to the observed 2-D velocity map from 
the FP data. As mentioned earlier, we adopt the 
kinematic projection angles from Barnes & Sell- 
wood (2003), which we also used for the photo- 
metric model. We therefore have to minimize x^ 
for the remaining free parameters, which are the 
mass-to-light ratios for the disk Tdj bulge Tb,i 
(if any), and parameters of the adopted halo. 

The line-of-sight velocity at the fcth-point in the 
projected disk plane predicted by our model is 

sill / cos / (•{)s(() — c)o) 

Vp,k =Vs+ Vc,T I (4) 



•^1 — sin^ i cos^((^ — (po) 



where Vc,t is the total circular speed from all 
mass components combined at that radius, and (j)o 
and Vs are, respectively, the previously-determined 
major axis position angle and systemic velocity. 

We evaluate the quantity Zk = {vd,k — Vp^k)/<^k 
at every pixel fc, where Vd.k is the estimated line- 
of-sight velocity derived previously by PW from 
the FP data. The quantity ak is the statistical 
uncertainty, estimated by PW from the Voigt pro- 
file fit to the FP data cube, plus a small additional 
constant term of 7 km/s added in quadrature to 
represent the likely velocity dispersion of gas in a 
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disk.^ The usual reduced Xr is therefore 
1 ^ 

fc=l 

where A'' is the number of pixels in the velocity 
map where we have an estimate of the line-of-sight 
velocity Vd.k, and the number of degrees pf free- 
dom, — N— the number of parameters. We 
omit only those pixels where the Voigt-fitter failed 
to find a statistically significant velocity estimate 
from the FP velocity scan. 

The value of Xr after minimization was unac- 
ccptably large for many galaxies. Inspection of the 
residual maps in these cases indicated that a small 
fraction of pixels, generally in the outer parts of 
the model, had differences of many between the 
model and data. The estimated velocities in these 
pixels generally differed significantly from those in 
the surrounding pixels, which is unlikely to be cor- 
rect on physical grounds, and such outlying points 
occasionally excrtcid undue influence on the fit- 
ted parameter values. The spurious velocities in 
these pixels generally arise from falsely identified 
"emission lines" that formally have small errors, 
which could be caused by incompletely removed 
cosmic rays, for example. Non-circular streaming 
motions, particularly in bars, arc further sources 
of major disagreement between our circular flow 
model predictions and the data. 

We therefore adopt a robust estimation proce- 
dure using Tukey's biweight, which reduces the in- 
fluence of non-normally distributed errors on the 
fitted parameters (Press et al. 1992). In this case, 
we minimize the function 

N ( „ »6 

^2 = ly-hl-^ + ^, \z,\<c 
k=i I T' otherwise 

where c is a constant. While data values that differ 
from the model prediction by < cakjl contribute 

^The intensity profile for any one pixel is broadened by the 
instrumental resolution, the intrinsic line- width from an in- 
dividual HII region, and by the different velocities of the 
various HII regions along the line of sight. The Voigt fit cor- 
rects for the instrumental broadening, which has a FWHM 
~ 125 km/s. An estimated velocity uncertainty < 7 km/s 
is therefore unreasonably low, and would give that datum 
too high a statistical weight in the fit. Furthermore, a sin- 
gle bright Hll region that dominates the spectrum in one 
pixel may not be as close to the circular speed as indicated 
by the fitted error if it were in the wings of the turbulent 
velocity spread of the ISM. 



to x& almost exactly as for the conventional Xrj 
the contribution from data values that are farther 
from the model prediction by ctTfe does not change 
as the parameters arc adjusted, and outlying data 
therefore do not influence the resulting parame- 
ter values. Press et al. (1992) recommend c = 6, 
but non-circular flow patterns in our galaxies sug- 
gest that a larger value would be more appropri- 
ate; after some experimentation we found c = 10 
still eliminated the influence of the extreme outly- 
ing velocities while allowing the great majority of 
the pixels to contribute to the fit with almost full 
weight. 

Since all values of Zk are adjusted during the 
minimization, the impact of the biweight changes 
at every iteration. It is therefore inadequate to 
first derive a 1-D rotation curve to which we can 
fit our models; we must always fit to the entire 
2-D velocity field. 

3.4. Steepness of the Inner Rise 

Many authors limit the upper value of T so 
that the model rotation curve rises no more steeply 
than the observed curve in the inner parts. We, on 
the other hand, try to minimize the difference be- 
tween the model and all the kinematic data with- 
out this constraint. As a result, our model rota- 
tion curve may exceed that observed in the central 
region. 

We do not regard mild overfitting in this re- 
gion as a serious shortcoming of our models. The 
reliability of the innermost data points has been 
discussed extensively in the context of halo con- 
straints from LSB galaxies (van den Bosch et al. 
2000; de Blok, McGaugh & Rubin 2001), and it 
is argued that every uncertainty leads to the ob- 
served velocity being lower than the true circular 
velocity. Non-circular motion, pressure support, 
seeing, and patchy emission can all lead to mild 
underestimation of the circular velocity in the ob- 
servational data from the inner galaxy. (We have 
no potential error from slit misplacement, because 
the position of the rotation center is evident from 
our 2-D velocity maps.) Tests to model worse see- 
ing gave results practically indistinguishable from 
those obtained with the original seeing. Patchy 
emission that is blurred by the seeing causes faint 
emission with the wrong velocity at points where 
no signal should have been detected. 
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Of these further sources of uncertainty, non- 
circular motions appear to have the greatest im- 
pact. Most galaxies in our sample have spiral pat- 
terns and some are clearly quite strongly barred. 
We have not, in this work, attempted to model 
the non-circular motions induced by such non- 
axisymmctric features, and continue to search for 
the best-fitting model with a circular flow pattern. 
This simplifying assumption clearly boosts xl > but 
it can also introduce systematic errors. In partic- 
ular, streaming motions in a bar generally cause 
the observed orbital speeds to be lower than circu- 
lar, because the gas spends more time, and is more 
easily detected, near the apocenter of the stream- 
ing pattern where the speed must necessarily be 
lower than circular. 

For these reasons, we consider a mild overfit 
in the central regions an acceptable result, and we 
accept a more significant overfit in galaxies having 
bars or other strong non-axisymmetric features. 

4. Stars-only Fits 

We first fit the velocity maps with models based 

on the light distribution alone (stars-only fits). 
The fitting parameter is Tjjj, with a second ^b,i 
parameter when a bulge is included. We list the 
best-fit parameters along with their uncertainties 
(estimated as described below) in Table 2. 

Figure 1 shows our results for four representa- 
tive galaxies. The mean orbital speed as a func- 
tion of radius is shown by the points with error 
bars; these values are deduced from the 2-D FP 
velocity map by the method described in Barnes 
& Sellwood (2003). The solid line in these Fig- 
ures shows the rotation curve of our best-fit, stars- 
only, model; the disk contribution is the dashed 
line, and the bulge contribution (if present) is the 
dash-dotted line. The four panels in Figure 1 show 
examples of different types of behavior. 

PW found that the large majority of galaxies in 
their sample could be modeled adequately with- 
out dark matter. Not surprisingly, we have also 
found that reasonable Tb,i and Td,i values can 
describe the data well in 34 of our 40 galaxies; 
Figure 1(a) shows a good example (ESO 268g44). 
In some cases {e.g., ESO 438gl5 and ESO 502g02) 
the model predictions begin to drop below the data 
in the outer few points, but match the shape of the 
inner data. 



The fit in the inner parts of eleven out of 
these 34 galaxies is significantly worse, due to 
bars and/or strong spirals that invalidate our as- 
sumption of circular orbits, but the amplitude and 
shape of the rotation curve over some radial range 
outside the barred region is predicted quite well 
by our mass model. Figure 1(b) shows an example 
of this behavior (ESO 439g20) - the model rota- 
tion curve runs significantly above the data values 
in the barred region, which is typical of these 11 
cases. 

Five of the remaining six cases exhibit the be- 
havior seen for ESO 444g47 (Figure Ic), in which 
there is no significant radial range where the pre- 
diction matches the data. The predicted rotation 
curve shape begins to drop well before the op- 
tical edge, whereas the data do not; this shape 
difference is a clear indication of a mass discrep- 
ancy within the visible disk. With no dark matter 
component available, the fitting routine returns an 
unreasonably large T/ in order that the underpre- 
diction farther out is balanced by a compensating 
overprediction of the observed speed in the inner 
parts. Since the T/ is larger than can be allowed 
by the inner rotation curve in these five cases, we 
follow PW and eliminate the outer data to find the 
maximum allowable T/. Subsidiary fits in these 
five cases are also reported at the end of Table 2. 

The remaining galaxy, ESO 381g05, shown in 
Figure 1 (d) , has large variations that are not well- 
predicted by the photometric model. This galaxy 
has a single strong spiral on one side, and may be 
disturbed by a recent, or on-going, merger. 

The total stellar masses we derive agree quite 
well with those found by PW from the same data. 
The agreement is not perfect for several reasons: 
In all but five cases, our fits include all the kine- 
matic data whereas PW limited their fits to the 
inner parts for all galaxies with outer mass dis- 
crepancies. Our disk-bulge decompositions are not 
identical and we also required Tb.i > O.ST/jj 
whereas PW frequently assigned much less mass 
to the bulge light. But the largest differences arise 
because we use kincmatically defined inclination 
angles, whereas they preferred the inclination de- 
rived from the photometric image; in a few cases, 
this difference exceeds 10°, leading to quite sub- 
stantial differences in Tj. In fact, the largest dis- 
crepancy is for ESO 381g05, which may be tidally 
disturbed, calling the assumption of an intrinsi- 
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cally flat disk into question. 

4.1. Uncertainties 

The uncertainty in the fitted T/ values is due 

mostly to the uncertainties in the galaxy's incli- 
nation, position angle, and distance. (Since we fit 
our models to a large number of data points, the 
T/ values are tightly constrained once these pa- 
rameters are fixed.) In a previous paper (Barnes 
& Sellwood 2003), we estimated systematic errors 
in inclination and position angles caused by spi- 
rals and other nonaxisymmetric structure in the 
galaxies. The best-fit value plus the estimated 
high and low limits for each of the two projection 
angles yield nine different combinations of inclina- 
tion and position angle which give nine separate 
estimates for T/ values. We take the uncertainty 
in an T/ value to be half the range of these nine 
Tj values. 

Since T is inversely proportional to distance, 
errors in our Hubble distances due to peculiar mo- 
tions will also bias our estimated values. Indeed, 
this sample of galaxies was originally observed to 
study peculiar velocities in the region toward the 
"Great Attractor" and may therefore have larger 
than usual peculiar velocities. Since typical red- 
shifts are 4000-5000 km/s, relative distance errors 
are unlikely to exceed 25%, however. Bothun et 
al. (1992) use the TuUy-Fisher relation to estimate 
distances, finding differences between redshift dis- 
tance and TF distance of typically about 25%. 
They find that the majority of these galaxies are 
actually closer than their Hubble distances, imply- 
ing that our T/ values may be underestimated by 
this factor. 

4.2. Discussion 

For our sample, the average Tdj is 2.6 and the 
average T b.i is 2.8 (using the T/ values of models 
with the best-fit i and (po). The average ^ d,i and 
f s,7 systematic uncertainties are std,/ = 0-5 and 
sr'nj = 0.7. 

Simulations by Bell &i de Jong (2001) predict 
that the mass-to-light ratio of a stellar population 
should correlate with its intrinsic B — R color. Bell 
et al. (2003) have applied the Bell & de Jong mod- 
els to a large sample of galaxies with multi-band 
photometry and find a somewhat shallower rela- 
tion, with a large scatter, which we show by the 



line in Figure 2(a). We plot our stars-only esti- 
mates of but it should be noted that these 
dynamical estimates are only upper limits - the 
stellar Toj could be lower if dark matter con- 
tributes significantly to the central attraction in 
the inner galaxy. The error bars on / are the 
systematic uncertainties discussed above, and the 
arrows indicate how ^ d,i would change if the TF 
distance were used in place of the Hubble distance. 

As we have only /-band images, we use galaxy 
colors obtained from the ESO-LV aperture magni- 
tudes, as recorded in NED. In principle, we should 
compare the T with the color of the disk alone, 
but NED gives magnitudes of the disk and bulge 
combined. The bulge light in our mostly late-type 
galaxies is a small fraction, typically < 10%, of 
the total in the J-band, so any correction for the 
bulge to the total color is likely to be small. We 
have corrected the NED magnitudes for Galactic 
extinction using values from Schlegel, Finkbeiner 
& Davis (1998) and for internal extinction using 
the prescriptions of TuUy et al. (1998). Errors in 
our estimated colors might be ^ 0.15 mag. 

Our average '^d,i values (and those of PW) are 
in reasonable agreement with those found by Bell 
et al. (2003) and it is interesting that there are no 
points significantly below their trend. A point ly- 
ing above the line, on the other hand, could simply 
indicate a significantly submaximal disk. While 
this could be the reason we see no trend, we first 
explore whether other effects could mask a trend 
in stars-only T/jj with color. 

It is extremely unlikely that the corrections 
for the bulge light to the total colors could be 
large enough to destroy a trend, and the galax- 
ies with no bulge (square symbols) also do not 
show one. Unfortunately, the number of galaxies 
without bulges is too small for a 2-D Kolmogorov- 
Smirnov test to determine whether this group 
is distributed differently from the bulge- and- disk 
galaxies in (T/, color) space. 

It is possible that a trend is masked by large 
errors in our values. The error bars include 

our estimate of the systematic error arising from 
uncertainties in the projection angles but do not 
include possible distance errors. They would have 
to be very large to mask the trend entirely, or con- 
spire to correlate with color, which seems imlikely. 

Bell et al. (2003) find the T values of blue stellar 
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populations are more sensitive than red to metal- 
licity. It is possible that the bluer galaxies in our 
sample are all metal-rich, but metallicity corre- 
lates with luminosity and we see no obvious trend 
of luminosity with color. 

We are suprised not to see any trend of stars- 
only Tu j with color, since a trend has been seen 
in other work {e.g., McGaugh 2004). The absence 
of a trend may imply that the halo contribution 
must correlate with color, which seems contrived 
for these bright HSB galaxies. The points in Fig- 
ure 2(b) are encoded in many ways to distinguish 
galaxies by luminosity, field galaxies from those in 
clusters, the membership of diS^erent clusters, and 
the arrows show where the point would move if we 
adopted photometric inclinations, instead of kine- 
matic. None of these appears to be a factor that 
could account for the absence of a trend. 

5. Dark Matter Halo Fits 

We next report fits to the velocity map that in- 
clude a dark matter halo component. Since our 
adopted halo models have two free parameters, 
we now have to minimize xl for cither three or 
four parameters, depending on whether a separate 
Tb,i is used. 

Estimating statistical uncertainties from the xl 
surface in a 4-parameter hypervolume is techni- 
cally more difficult than for the 2-parameter stars- 
only models. The difficulty is compounded be- 
cause we often find that the adopted halo pa- 
rameters are strongly degenerate - that is, the 
minimum in the Xb hypcrsurfacc lies in a long, 
curved, and almost flat-bottomed valley. We es- 
timate the uncertainties in each fitted parameter 
using a Markov Chain Monte Carlo method {e.g., 
Christensen et al. 2001; Kosowsky, Milosavljevic & 
Jimenez 2002), deriving the estimated likelihood 
of a set of parameters from the xl value. Our 
Markov Chains for each galaxy contain 1 million 
elements and appear to be well-converged. We 
construct a histogram of the chain elements in 
the multi-dimensional parameter space; the like- 
lihood region is the set of all bins with more than 
some threshold number of chain elements, with the 
threshold chosen such that the selected bins con- 
tain 68% of the chain elements. The set of chain 
elements in the selected likelihood region can then 
be projected onto the parameter axes to determine 



the statistical uncertainties. Figure 3 shows an ex- 
ample for the Pseudo-isothermal (PI) halo model 
of ESO 322g82. Strictly speaking, these confidence 
regions are not exactly the 1-a imcertainties, since 
we define likelihoods from x^, and not the usual Xr 
function. However, we find in practice that work- 
ing with Xrj cases where the outlying points 
have little effect on the position of the minimum, 
leads to very similar confidence regions. 

We use an F-test to determine whether a halo 
is required to fit our data. This test determines 
whether the reduction in xl min from the stars only 
case to when the fit includes the extra two halo pa- 
rameters is large enough to justify the additional 
parameters. (Again, we treat xl ^ adequate 
surrogate for Xr ^'^d, of course, we always com- 
pare with Xr the entire velocity map from the 
stars-only fits.) We discard the fits in a few cases 
only where this test indicates that there is less 
than 95% chance that the halo is necessary, since 
the values of the parameters in such cases are un- 
likely to be meaningful. Just three galaxies (ESO 
317g41, ESO 382g06 & ESO 438g08) fail this test 
for all the types of halo we have fitted. 

We also find in many cases that the halo drives 
Tz3,/ to very low values, and the observed flow pat- 
tern is determined largely by the halo. It is clear 
that Td,i <C 1 is unphysical - the stellar popula- 
tion of the disk clearly must have some significant 
mass. Since the galaxies in our sample have quite 
normal colors, we regard Tdj < 0.5 to be unre- 
alistically low, and therefore discard any fits for 
which the 1-a upper bound on Tjjj < 0.5. Very 
few cases fall just below this cut, and T oj <C 0.1 
for the large majority discarded for this reason. 
This criterion leads us to discard one galaxy in 
our sample (ESO 323g39) for all types of halo we 
have tried. 

We have not attached any significance to the ac- 
tual Xb min value, which always exceeds unity and 
is mostly greater than 2 ranging up to greater 
than 11 in one case! Formally at least, large val- 
ues of this parameter indicate that the model is an 
unacceptable fit to the data, which we do not dis- 
pute. However, examination of the residuals in the 
2-D velocity map indicates that the large xlmm 
values arise from three principal sources: outlying 
points, bars, and spirals, in decreasing order of 
importance. We argue that large values of xl min 
do not invalidate our estimates of halo parameters. 
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The biwcight function is a robust statistic designed 
to minimize the influence of the outlying points, 
even though they still make a large contribution 
to the value of Xb min i'^'i- ^) ■ ^^"^ h.alo parame- 
ters are determined more by the overall shape of 
the rotation curve than by the non-axisymmetric 
streaming motions, which could possibly be mod- 
eled with a lot of effort {e.g. Weiner et al. 2001b), 
but do not affect the final parameter values to a 
great extent (Weiner et al. 2001a). 

5.1. Pseudo-isothermal Halos 

The PI density distribution is 
Po 



Ppi = 



l + (r/re)2' 



(7) 



specified by two parameters po, the central density, 
and Tc, the core radius. The associated squared 
circular velocity is 



^c,PI 



1 arctan — 



(8) 



where v^c = (47rG/9o)^^^'"c is the asymptotic circu- 
lar velocity. 

The best-fit parameters and their associated 
uncertainties for all 40 galaxies are given in Ta- 
ble 3. In only 3 cases out of the 40 galaxies in our 
sample does the addition of the extra two param- 
eters for the halo lead to such a small reduction 
in Xftmin that the halo is not needed; the value 
of Xb min these cases is marked with an aster- 
isk; we discard these three cases in the following 
discussion. 

Ten of the remaining 37 galaxies have T u j val- 
ues for which the 1-a range never exceeds 0.5. 
Non-circular motions in the bar may be one of 
the factors that drives the solution towards a low 
Tj value and a dominant halo. In these cases, the 
circular speed is often underestimated, while the 
rotation curve predicted by the luminous matter 
has a more steeply rising shape. Thus the fit can 
be improved by decreasing the T/ and fitting the 
data mostly with the halo, which has a more ad- 
justable shape. Whether this happens depends on 
other factors also, such as the bar orientation. 

The halo parameters {po and Tc) of the remain- 
ing 27 galaxies span a wide range, as shown in Fig- 
ure 4. The line of slope —2 in this plot is not a fit 
to the data, but shows the relation = AnGporl 



for Woo = 200 km/s. Galaxies falling above (be- 
low) this line have larger (smaller) values of t^oo- 
The fitted points are not randomly distributed 
about the line: the points above (below) the line 
tend to have large (small) Tc values. The distribu- 
tion suggests, perhaps only marginally, that larger 
Voo galaxies, which also tend to be more lumi- 
nous, have larger, lower-density cores than do the 
smaller galaxies. Despite using a slightly differ- 
ent cored dark matter distribution, the 6 galax- 
ies presented in Gentile et al. (2004) follow es- 
sentially the same relation between central den- 
sity and core radius. Donato, Gentile, & Salucci 
(2004) also present isothermal fits for 25 galaxies. 
Ten are "regular" spirals (neither dwarf nor low 
surface brightness), and they too lie among our 
points in the (poj^'c) plane. Donato, Gentile, & 
Salucci (2004) have additionally found a tight cor- 
relation between the disk scale length and Tc. PW 
fit exponential disks to their galaxies and deter- 
mined disk scale lengths {Rd), and we adopt these 
values as reasonable estimates. Despite the trend 
for the galaxies with larger v^, and presumably 
larger disks, to have larger Vc values, we do not re- 
produce their claimed correlation between i?^ and 
Tg. Kormendy & Freeman (2003) find pq oc rj^'^ 
for PI halo fits, which is a shallower slope than 
the trend shown our Figure 4. However, their fits 
are made with the assumption of a maximum disk, 
which precludes central densities as high as some 
that emerge from our fits. 

The three error bars that extend from points 
well beyond the right-hand edge are from three 
cases (ESQ 322g45, ESQ 323g73 & ESQ 375g02) 
in which halo core radius is huge and very poorly 
constrained; the halo contribution to the circular 
speed rises essentially linearly to the outermost 
observed velocity. We are therefore unable to de- 
termine Voo with any confidence. Such degener- 
acy has been reported by other authors, even by 
Broeils (1992) who had high quality HI data ex- 
tending to radii well outside the optical disk. 

5.2. NFW Halos 

The collapse of dark matter halos in simula- 
tions of large-scale structure formation yields halo 
density profiles that are not at all like the pseudo- 
isothermal function used above. The precise form 
of the density profile that emerges in these simu- 
lations is still disputed (Navarro, Frenk & White 
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1996; Power et al. 2003; Dicmand, Moore & Stadcl 
2004), but all agree that the density increases 
steeply towards the center. We have therefore 
tried a broken power-law for the halo density pro- 
file. In this work, we treat the halo as spherical 
and make no attempt to take account of halo com- 
pression as a result of the cooling and settling of 
baryons to form the luminous component. Our 
objective is simply to assess whether the data are 
more consistent with a cusped density profile in 
galaxies today, and not to test detailed predictions 
from cold dark matter (CDM) simulations. 

We adopt the original NFW density profile: 



Pnfw 



Ps 



r{r + TsY 



(9) 



where ps is a scale density and is a scale length. 
The squared circular velocity arising from this 
density profile is 



^c,NFW 



= V. 



200 



ln(l -Hx) -x/{l+x) 



ln(l + c) - c/(l + c) 

(10) 

where x = r/rg, c = r^Qo/rg is the concentra- 
tion parameter and W200 = lOcrg/fo is the circular 
speed at r2oo- Inside the radius r2oo the average 
density is 200 times the critical density of the uni- 
verse (3if^/87rG). 

The best- fit NFW parameters (c and v^m) are 



listed in Table 4. The reduction in 



is too 



small to justify the two halo parameters in 4 cases 
out of the 40 galaxies in our sample; the value of 
Xftmin these cases is marked with an asterisk. 
Once again, we discard these for the remainder of 
this section. We also discard 16 of the remaining 
36 galaxies because the 1-a range of their Td,i 
values never exceeds 0.5. 

As with the pseudo-isothermal halo fits, the 
best fit NFW halo parameters of the remaining 
20 galaxies have a considerable range, as shown 
in Figure 5. In only 9 cases, which tend to have 
higher values of c, the 1-ct range of ^^200 does not 
extend up to values greatly in excess of any ob- 
served velocity. The enormous cirror bars for the 
other 11 galaxies indicate that the halo parameters 
are strongly degenerate. The line is the mean cor- 
relation between these halo parameters in a con- 
cordance ACDM simulation (Bullock et al. 2001). 
While our model halos are not directly compara- 
ble to those from simulations (see § 6.3 for a more 



thorough discussion of the comparability of the 
simulated and model halos) and the uncertainties 
are quite large, the points in Figure 5 appear to 
describe a steeper trend than the prediction. 

5.3. Power- law Halos 

The results from fitting the above two standard 
halo models were not particularly informative for 
a number of reasons. Extreme minimum disks are 
often favored in a significant fraction of the galax- 
ies (mostly those with strong non-axisymmetric 
features), which we then discard as unphysical. 
Also, the halo parameters are frequently highly 
degenerate since the minimum of the hypersur- 
face lies in a long, narrow, curved valley. 

However, rotation curve data do tightly con- 
strain something, namely the mass or mean den- 
sity interior to the last measured point. The just- 
noted degeneracy in the NFW parameters arises 
because halos having the allowed wide range of pa- 
rameters all have closely similar density profiles in 
the inner region that is the only part constrained 
by our data. Cases where no mass is assigned to 
the luminous material appear to arise because fit- 
ting functions with a length scale and a density 
scale as free parameters have too much freedom 
to adjust to quirks in our data. 

We have therefore also tried a halo density pro- 
file with a smaller range of possible shapes: a sin- 
gle power law with no characteristic length scale. 
The density profile we adopt is 



Pa(r-)=P((^) , 



(11) 



where pi is the halo density at r;, the radius of 

our outermost velocity measurement. We stress 
that we fit this function only to the radial range 
spanned by our data and cannot say how the den- 
sity profile continues beyond this radius. The 
squared circular velocity in this halo is 



2 2 



n 



(12) 



Again we have two free parameters, Vi and a, with 
pi = {Z-a)vy{A-KGr}). 

Table 5 gives the best-fit parameters and uncer- 
tainties for the power law fits. The reduction in 
Xfc mill is too small to justify the two halo parame- 
ters in 6 out of the 40 galaxies in our sample; the 
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value of xl in these cases is marked with an as- 
terisk. Once again, we discard these. We need to 
discard only 2 of the remaining 34 galaxies because 
the 1-(T range of their Toj values never exceeds 
0.5. This small fraction contrasts with the 10 and 
16 cases for the PI and NFW fits respectively. 

With this halo form, we do not find strong de- 
generacies between the fitting parameters, imply- 
ing that the fitted values are much more robust. 
Figure 6 shows that the halo parameters are also 
spread over a more limited range. Yet despite the 
more limited range of possible shapes, the values of 
Xfe min ^'''^ scarcely significantly larger than for the 
more flexible functions, and are actually smaller 
in a few cases. 

The reason the power-law halo model is more 
useful, as judged by the preceding statistics, is 
that it has less flexibility. The length scale in 
the pseudo-isothermal and NFW halo models can 
produce a feature in the halo contribution that 
matches a feature in the data, thereby yielding 
a lower Xr^ even though this is often achieved at 
the expense of driving Tdj to an absurdly low 
value. In these cases, it is likely that some non- 
axisymmetric disk feature drives the fit to this part 
of parameter space, rather than a true feature in 
the halo profile. 

Figure 6 shows that the distribution of a val- 
ues is concentrated between 0.8 < a < 1.6, al- 
though a few values are close to zero and some 
range up to a < 2. The significance of a > 1 
is unclear, however, since the probable halo con- 
tribution should be something in between uniform 
rotation and exactly constant orbital speed, which 
requires < a < 2. While all intermediate values 
yield a halo rotation curve that rises fast initially 
and bends over to a less steeply rising part towards 
the outer radius, the abruptness of this transition 
does increase with a, and the concentration of val- 
ues in the range 0.8 < a < 1.6 may be significant. 
Recent cosmological N-body simulations have also 
found similar a values (Diemand, Moore & Stadel 
2004). 

6. Disk Masses and Halo Densities 
6.1. Mass-to-light Ratios 

We plot the Tdj values of the three halo mod- 
els versus color in Figure 7. The error bars now 
reflect only the statistical uncertainties in Tn,! 



values, and do not include any estimates of depro- 
jection uncertainties which are of comparable mag- 
nitude, but which would be burdensome to com- 
pute. Since the colors are the same, the abscissae 
and their uncertainties in each panel are the same 
as in Figure 2. Each panel of this figure has fewer 
than 40 points because we have discarded any fit 
for which the halo was unnecessary or that yielded 
such a low Tdj value that its 1-a range never ex- 
ceeded 0.5. We again indicate the relation fitted 
by Bell et al. (2003). The most striking feature of 
these plots is that adding halos has not changed 
the flat character of the distribution of points seen 
in Figure 2; we still do not see a trend with color 
for any of the three halo forms. 

Our velocity data in a few cases do not re- 
quire a halo component at all. This statement 
is true for all halo types in just three galaxies: 
ESO 317g41, ESO 382g06, and ESO 438g08. ESO 
317g41 and ESO 438g08 have strong spirals, but 
ESO 382g06 is axisymmetric. The kinematic data 
for ESO 317g41 extends only to 80% of R23.5, the 
radius at which the /-band surface brightness is 
23.5 mag/arcsec^. This raises the possibility that 
the radial extent of our data too small to probe the 
halo. However, for ESO 382g06 and ESO 438g08, 
the kinematic data extend at least to R23.5, ap- 
proximately 4 disk scale lengths. The rotation 
curves of these two galaxies turnover and decrease 
in the outer parts. As these galaxies appear to be 
baryon dominated, one might expect good agree- 
ment between our values and those predicted 
by Bell et al. (2003). In fact, we find that only 
ESO 438g08 is marginally consistent with the pre- 
dicted value while the other two are several sigma 
above the prediction. Even though none of the 
simple halo functions tried here improves the fits 
by enough to justify the extra parameters, these 
three galaxies may yet have dark matters halos. In 
fact; if the halo rotation curve has a similar shape 
to that of the disk in the range covered by our 
data, any relative fraction would yield an equally 
acceptable fit. 

Some galaxies, such as ESO 216g20 - the red- 
dest in our sample, do seem to require a halo that 
makes their fitted 7 fall below that expected 
for a stellar population of its color. We note that 
the TF-distance for this galaxy is larger than its 
Hubble distance, so its unusually low is un- 

likely to result from too short an adopted distance. 
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ESO 322g44 is another case which falls below the 
line in all three panels and there are others which 
are low when we accept the fit, and are discarded 
as unrealistically low for other halo types. 

Of the galaxies lying above the line from Bell 
et al. (2003) in Figure 2, only one of the bluest 
galaxies in our sample. ESO 446g01, remains sig- 
nificantly above the lino in all the T/ vs. color 
plots, while three other galaxies (ESO 381g05, 
ESO 437g31, ESO 501g01) lie well above the hne 
for the type(s) of halo(s) for which we accepted 
the fit, but are discarded for other types. The in- 
troduction of any of three halo types into these 
four galaxies did not reduce their Td,i values to 
the predicted value. 

None of our adopted halo models yields Td,i 
values that follow the trend found by Bell et al. 
(2003). It is possible that systematic errors such 
as distance and/or inclination errors, radial ex- 
tinction variations, gas mass corrections, thick- 
ness corrections, non-circular motions, etc., have 
masked the expected trend, but they would have 
to be as large as factors of three, and in a sense 
to correlate with color. It is more likely that our 
models have mis-estimated the dark matter con- 
tent in the majority of these galaxies, and our 
'^D,i values are spurious. 

Figure 8 compares the '^d,i values obtained 
from each of our fits to the 2-D velocity maps. 
Each panel displays considerable scatter. The val- 
ues from adopting one halo model are frequently 
inconsistent with those from another model at the 
several-a level, although a few cases, which gen- 
erally have smaller uncertainties, do lie nicely on 
the diagonals in all three panels. The open circles 
mark galaxies which fail the F-test for either of 
the halo models. The most significant feature in 
this Figure is that the NFW halo density profile 
generally leads to the lowest assigned '^d,! value 
of the three functions we have tried. 

6.2. Halo Masses 

We have determined the total mass in the halo 
interior to our outermost available velocity mea- 
surement. The halo masses for the three differ- 
ent functional forms, normalized by total enclosed 
mass, are compared in Figure 9. Again, the open 
circles mark galaxies which fail the i^-test for ei- 
ther of the halo models. We have used our MCMC 



study to determine the confidence intervals, which 
range up to 70% in the worst few cases but are 
generally much lower. 

It is clear from this Figure that the NFW form 
tends to require a larger halo mass fraction than 
the other two types, which is consistent with its 
generally lower disk masses. Many cases require no 
halo or negligible disk mass, and other cases show 
large differences in fitted halo masses between the 
different halo models. 

In general, the agreement between the halo 
masses calculated for the different halo types is 
poor. For 12 galaxies, the total uncertainty in 
PI halo mass fraction is < 0.2. These galaxies 
are, on average, 2.3(7 away from agreeing with the 
NFW halo mass fractions and are 2.1a away from 
agreeing with the power-law halo mass fractions. 
The 5 galaxies which have NFW halo mass frac- 
tion uncertainties < 0.2 are l.Scr and 2.1cr from 
agreeing with the PI and power-law values, respec- 
tively. Finally, 17 galaxies have power-law halo 
mass fraction uncertainties < 0.2. These galax- 
ies are 2.2(T and 2.3cr away from agreeing with the 
PI and NFW values, respectively. However, in a 
few cases the halo mass fraction has a small un- 
certainty, is neither nor 1, and lies close to the 
diagonal in this Fig. 9. One such case is ESO 
445gl5, for which we claim a consistent and well- 
constrained halo mass fraction for all three density 
profiles we tried. 

Aside from these rare cases, we conclude that 
not even our 2-D data constrain the dark mat- 
ter content of a galaxy - the fitted amount is ei- 
ther highly uncertain, or depends on the halo mass 
function adopted. 

6.3. Halo Densities 

Cosmological simulations make predictions for 

the central densities of dark matter halos. Specifi- 
cally, Alam, Bullock & Weinberg (2002) introduce 
a measure of the central densities of halos, A„/2) 
which is the average halo density inside the ra- 
dius at which the velocity is half the maximum 
halo velocity [p{Ry^i2)], normalized by the critical 
density of the imiverse, pc = 3HQ/{8nG). Their 
Figure 2 plots the predicted from an ACDM 
simulation versus maximum halo velocity Wmax in 
which larger dark matter halos (high Wmax) have 
lower densities. The same figure has points de- 
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rived from real galaxies which do not appear to 
follow the predicted trend. 

We show our derived A„/2 versus Vmax for the 
NFW halos only in Figure 10. The line illustrates 
the trend expected from simulations (Zentner & 
Bullock 2002). One point marked by an open cir- 
cle is one of the four cases for which the F-test 
indicates no halo is required, two other cases are 
points off the left side of the plot, and the fourth 
case cannot be plotted as the best fit assigns zero 
halo mass. The points without symbols, which lie 
mostly above the line, are fits with very low Td,i- 
The rest are marked with filled symbols. Those 
that have enormous error bars are cases that are 
degenerate, because our optical data do not extend 
far enough to constrain r^. 

Our data points are not directly comparable 
with predictions of ACDM, since the simulated 
halos should be compressed by baryonic infall. To 
make a detailed test, we should include halo com- 
pression as a step in the minimization loop in or- 
der to determine the parameters of the uncom- 
pressed halos, possibly also testing a range of un- 
compressed density profiles as part of the uncer- 
tainty in the predictions. Our less ambitious ap- 
proach is simply to test whether the NFW profile 
describes galaxy halos in the local universe any 
more closely than the others. Thus our parame- 
ters, if they have any relevance at all, are for com- 
pressed halos. It is clear that the mean densities 
of the uncompressed halos would be lower, by an 
amount that depends on the stellar mass. 

The large error bars on many points make it 
difficult to conclude much from this Figure; we 
cannot claim that our data are inconsistent with 
the prediction, but neither do they offer much sup- 
port for it. This result is similar to the outcomes 
in other studies that compare simulations to pa- 
rameters derived from galaxies (Alam, Bullock & 
Weinberg 2002; Zentner & Bullock 2002). Even 
points with small error bars have a large scatter 
about the line, whereas the simulations predict a 
rms scatter of about a factor 3 (Alam, Bullock 
& Weinberg 2002). We also note that as points 
without symbols represent fits that attribute es- 
sentially no mass to the luminous matter, more 
realistic fits would attribute a lower average halo 
density and move those points to lower A^,/2 val- 
ues. 

We cannot make a similar plot for the power 



law halo, since we know nothing about how the 
profile may vary at radii beyond our data. How- 
ever, we did obtain the lowest halo mass with this 
fmiction (Fig. 9) which, combined with the greater 
compression from the more massive disks suggests 
that the densities are considerably lower than for 
the NFW halos. 

7. Rotation Curve Fitting 

In this section we compare our results from the 
2-D velocity maps with what could be achieved 
for the same galaxies from 1-D longslit data. As 
we do not have separate longslit data, we extract 
pseudo-data from the 2-D velocity map along a 
line across the major axis of the galaxy. 

7.1. Pseudo-slit Rotation Curves 

We mimic data that could result from longslit 
spectroscopy, by placing a "slit" along the pho- 
tometric major axis of each galaxy and passing 
through the point of brightest continuum emis- 
sion. We extract the mean velocity, vios, and its 
uncertainty from our 2-D velocity maps every 1.5", 
a typical seeing length, from the pixel whose cen- 
ter is closest to the chosen line; interpolation or 
averaging over several pixels is unnecessary, since 
the pixel sizes are generally less than the seeing. 
The S/N in any one pixel of our map is probably 
less than would be obtained from a longslit obser- 
vation, but not by a large factor. (The FP scan 
used 10 min exposures at each filter setting, but as 
the device used has about twice the throughput of 
a typical spectrograph, our S /N may be equivalent 
to that from a 20 min spectrographic observation.) 

We treat these data by the typical procedures 
for longslit data {e.g., de Blok, McGaugh & Ru- 
bin 2001). We estimate the systemic velocity by 
folding the fios curve, by eye, assuming that the 
receding and approaching sides are symmetric. 

Since a velocity map allows us to determine Vg 
with greater precision, we have investigated the 
differences between the Vg^map and Vg^snt values. 
The rms difference between the systemic veloci- 
ties derived from the pseudo-slit data and our ve- 
locity maps is Agyg = 10 km/s. It is interest- 
ing to break this comparison down for axisym- 
metric, spiral, and barred galaxies. Axisymmet- 
ric galaxies (disks without strong spiral structure) 
have Agys = 6 km/s. Galaxies with spirals have 
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Agys = 13 km/s. Overall, barred galaxies have 
Asys = 10 km/s. We further subdivide the barred 
galaxies according to the position angle of the bar. 
When the bar lies along the minor axis, Agys = 8 
km/s. Galaxies with bars at an oblique angle to 
the major axis have Agyg = 12 km/s and those 
with major axis bars have Agyg = 10 km/s. In the 
central regions of these galaxies where typically 
the rotation curves rise approximately linearly, the 
uncertainty in Ws.siit relates directly to an uncer- 
tainty in the position adopted for the center. Us- 
ing the maximum uncertainty quoted above leads 
to an uncertainty in dynamical center positions 
< 1". 

We retain all the points (receding and ap- 
proaching) as independent data points in the rota- 
tion curves. Finally, we correct the rotation speeds 
for inclination using the photometric value since 
this is all slit observers would have at their dis- 
posal. 

7.2. Comparison of Map and Pseudo-slit 

Rotation Curves 

Our pseudo-slit rotation curves generally agree 
well with those we derived from the 2-D veloc- 
ity maps. A few illustrative cases arc shown in 
Figure 11; the dotted error bars denote pseudo- 
slit data, while the solid points are from our fits 
to the 2-D maps. Panel (a) shows an example 
(ESO 322g82) of good agreement between the two 
curves. The case in panel (b) shows a galaxy with 
strong spirals (ESO 322g36); velocities differ be- 
tween the approaching and receding sides over the 
radial range from 3 to 8 kpc, but the velocity 
map rotation curve nicely bisects this asymmetry. 
The agreement between the two curves in panel 
(c) (ESO 438g08) is poor, which results from a 
disagreement between the photometric major axis 
that a slit observer would use and the major axis 
that we find from our 2-D velocity map; the spi- 
rals present in this galaxy have biased the esti- 
mate of photometric position angle (see Barnes & 
Sellwood (2003) for a fuller discussion of this ef- 
fect). The case shown in panel (d), ESO 501g01, 
shows the greatest advantage of velocity map data 
over slit spectroscopy, since this galaxy has major 
axis Ha emission over only a limited radial range, 
whereas the velocity map allows a much greater 
radial range to be explored. 



7.3. Mass Models 

We have used our pscnido-slit rotation curves 
to determine T/ values in a manner similar to 
the 2-D data. We fit constant T/ models, to 
the inclination-corrected velocities, using the same 
procedure as described above, but for this much 
smaller number of pixels. As before, we first fit 
models without halos and then include the three 
different halo types. 

7.3.1. Stars-only Fits 

The comparison of pseudo-slit and map T/ val- 
ues is illustrated in Figure 12. The pseudo-slit and 
map values have rms fractional diff'erences of 
22% while the Tbj values differ by 39%. We find 
that in 28 of 40 galaxies the pseudo-slit val- 
ues are 25% lower than the map values. The 
remaining 12 galaxies have map T u j values which 
are 13% less than the pseudo-slit values. In 17 of 
the 31 galaxies with bulges, the pseudo-slit Tb,i 
values are 30% lower than the map Tbj values. 
The other 14 galaxies have map Tb,i values 49% 
smaller than the pseudo-slit values. A few outliers 
inflate these averages, so the overall Td,i differ- 
ences are closer to 15% and the Tb,i differences 
are more like 25%. Three of the worst outliers in 
the Tdj comparison (ESO 381g05, ESO 382g06, 
ESO 445g39) have quite poor agreement between 
the pseudo-slit and velocity map rotation curves. 
Two others (ESO 323g39 and ESO 438g08) have 
slit misalignment issues. However, the trend for 
pseudo-slit Tn,i values to be lower than the map 
Td,/ values remains. 

This trend may have an explanation in the dif- 
fering methods. Fitting pseudo-slit data gives 
equal weight to each data point regardless of ra- 
dial location. However, velocity map fits utilize 
more data as radius increases, giving more weight 
to the outer regions of the rotation curves. Since 
the outermost points tend to have higher velocity, 
the map values refiect this by being slightly 
larger. We find a similar bias when halos are added 
- see below. 

As expected for fewer data points, the average 
pseudo-slit statistical uncertainties are « 3 times 
larger than the velocity map statistical uncertain- 
ties. The systematic uncertainties due to position 
angle uncertainties arc of the same magnitude for 
both methods, and the estimated systematic error 
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in Tj w 0.5. As with the stars-only map fits, the 
T/ values do not show any trend with color. 

7.3.2. Halo Fits 

As with the 2-D velocity maps, we have fitted 
the three halo models to these rotation curves from 

1- D pseudo-slit data. Figure 13 shows how the 
halo mass fraction (within the radius of the outer- 
most fitted point from the velocity map) compares 
between the fits to the pseudo-slit data and the ve- 
locity maps. 

Not surprisingly, we find constraints on the halo 
masses are generally weaker for slit data than for 

2- D maps, by a factor > 2 in most cases. The sig- 
nificant outliers in these plots are generally cases 
of disagreement between the photometric inclina- 
tion adopted for the pseudo-slit, and the kinematic 
inclination used for the velocity map. 

The slit-estimated halo masses tend to be 
higher than the map-estimated values, which we 
attribute to the above-noted significance of the 
large number of pixels in the outer parts of the 
map. Since slit data give greater relative weight 
to the inner points, an overfit in the inner parts 
carries a greater penalty, sometimes forcing Tjyj 
to lower values, which requires more halo mass to 
match the outer data. 

8. Summary &: Conclusions 

Palunas & Williams (2000) present a sample 

of homogeneously collected, 2-D optical velocity 
maps and /-band surface brightness distributions 
of galaxies. We have re-analyzed a subsample 
of these galaxies, omitting mostly those that are 
highly inclined, in order to determine whether 
such data constrain the dark matter content of 
these galaxies any more tightly than do previ- 
ous studies using longslit data or velocity maps 
at lower spatial resolution. 

Our results echo those of PW in finding that 
the stars-only model fits the shape of the observed 
rotation curve over a significant radial range in 
34 out of the 40 cases. This fraction includes 11 
cases where bars or spirals cause a mismatch in 
the inner barred region which we argue is accept- 
able because of non-circular motion. A further 5 
cases could be fitted tolerably well when the outer 
data were discarded, as found by PW. In addition 
to these 5 galaxies, several others have noticeable 



mass discrepancies in the outer parts. These stars- 
only fits yield an average Tdj of 2.6 and an aver- 
age Ts J of 2.8, similar to the values found by PW. 
The main sources of uncertainty in these estimates 
are those arising from distance errors and possible 
errors in the adopted projection geometry. 

These average T/ values, even with no dark 
matter, are in the range predicted from stellar 
population models, but we do not find a trend of 
T/ with color. Instead of finding lower T/ values 
for bluer galax;ies, we find a fiat distribution which 
does not decrease for the bluest galaxies. 

We fitted three dark matter halo profiles; 
two are standard, pseudo-isothermal (PI) and 
Navarro- Frenk- White (NEW), and we have also 
used a simple power law function over the radial 
range of our data. Three galaxies in our sample do 
not require a halo of any type we tried, although 
these galaxies may have dark matter halos of some 
other form. Should the true inner halo rotation 
curve shape resemble that of the luminous matter, 
any relative fraction would yield an equally accept- 
able fit; studies such as this, which are based on 
rotation curve data alone, could never break such 
a degeneracy. 

In most cases, the fit to the data is signifi- 
cantly improved by adding one or more of our 
adopted halos, but we do not find a preference for 
any particular halo density profile, in the sense 
that the values of Xbmin rather similar 

no matter which form is adopted. However, we 
find wide variations in the T^j values and halo 
masses, within the radial range spanned by our 
data, depending upon which halo density profile 
is adopted. This finding implies that even with 
high-quality 2-D data, rotation curve fitting is un- 
able to determine either the disk mass-to-light ra- 
tio, halo mass, or density profile with any degree 
of confidence. 

The pseudo-isothermal and NFW halo models, 
which have both a length scale and a density nor- 
malization as free parameters, often lead to fitted 
Tj values that are absurdly low, which seems to 
indicate that such functions have too much free- 
dom to adjust their shapes to quirks in the velocity 
maps, such as may arise from non-axisymmetric 
features in the disk. Tlicise functional forms also 
lead to strong degeneracies in the fitted parame- 
ters, especially for low halo densities. Our sug- 
gested power-law halo density profile is more ro- 
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bust, in the sense that it does not have strong de- 
generacies and requires fewer near-massless stellar 
populations. 

Following the suggestion of Alam, Bullock & 
Weinberg (2002), we have calculated the central 
density parameter A„/2 for the fitted NFW halos. 
We find that our fitted halo densities are highly 
uncertain, but seem to span a very wide range 
that includes values consistent with the predic- 
tions from c;osniological simulations but also ex- 
tend to much lower densities; a few galaxies have 
no detectable dark matter at all. 

Fits with the NFW halo function generally at- 
tribute the least mass to the disk, and the most 
to the halo, while the power law prefers the high- 
est disk masses. On average, the power-law halos 
have radial density profiles oc or somewhat 
steeper, i.e. a slope resembling the inner parts of 
halo density profiles from recent cosmological sim- 
ulations of structure formation (Dicmand, Moore 
& Stadel 2004). However, the corresponding mean 
halo densities are generally lower than those ob- 
tained with the NFW function. 

We have compared our fitted Tjjj values with 
population synthesis models. While our stars-only 
T oj values (Fig. 2) arc in the range predicted by 
Bell et al. (2003), we do not find a trend with 
color. When halos are included (Fig. 7), many 
T oj values are lower than their predictions. Dis- 
tance errors affect the estimated T/ directly, but 
it is unlikely they could be large enough to bring 
the points into agreement with the prediction; fre- 
quently, adopting TuUy-Fisher distances instead of 
Hubble distances would reduce Tn,i still further. 
Also, wc do not sec a correlation between color 
and distance, luminosity (metallicity), or cluster 
membership needed to recover the predicted trend. 
The most likely explanation is that none of our 
models has captured the actual distribution of 
dark matter in our sample of galaxies and, with- 
out this knowlegde, our fits cannot determine the 
true Tn,!- 

We have compared our results from fitting 2- 
D velocity maps with those from pseudo-slit data 
extracted from our data on the photometric ma- 
jor axis only. Derived rotation curves generally 
agree reasonably well with the rotation curves 
from the velocity maps. In a few cases, dis- 
agreement stems from misaligned slits, poor ma- 
jor axis data, and/or nonaxisymmetric structure. 



Not surprisingly, we find the uncertainties in the 
estimated parameters are « 3 times larger from 
pseudo-slit data than from the 2-D maps. We find 
a tendency for T/ values from our pseudo-slit data 
to be f« 15% lower than the corresponding velocity 
map values, which we attribute to the extra weight 
that 2-D maps give to the larger number of data 
values in the outer parts of the map. The great- 
est improvements afforded by 2-D velocity maps is 
that they offer tighter constraints on the projec- 
tion geometry and the systemic velocity, as well as 
revealing more extensive rotation curves in cases 
with sparse data along their major axes. 

Fits of axisymmetric models to two-dimensional 
velocity maps, even at optical resolution, do not 
provide significant constraints on either the mass- 
to-light of the luminous matter or the halo mass 
fraction or its density profile. Additional data ex- 
tending to larger radii, e.g. from neutral hydro- 
gen maps, must tighten the constraints somewhat, 
but the improvement is not dramatic; e.g. Blais- 
Ouellette, Amram & Carignan (2001) still find 
large variations in Td,i and halo mass between 
the different halo models. We therefore conclude 
that additional information is needed to constrain 
the disk mass, and therefore also the halo density, 
in a galaxy. 
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Fig. 1. — Various types of rotation curves for the stars-only models. The dashed line represents the disk 
contribution. The bulge contribution is shown with the dash-dotted line. The solid line is the combination 
disk-|-bulgc. a) ESO 268g44: an example of good agreement between model and data, b) ESO 439g20: 
a barred galaxy for which the model overpredicts the inner rotation curve; the bar has a half-length of 
approximately 2 kpc. c) ESO 444g47: an example of an overfitted disk driven by an evident mass discrepancy 
in the outer regions, d) ESO 381g05: a galaxy with strong spirals and a possible tidal distortion. 
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Fig. 2. — Fitted Td,i values for stars-only vs. dereddened B — R. In both panels the line is the best-fit 
correlation found by Bell et al. (2003). (a) Estimated systematic errors arc shown and the arrows indicate 
how the points would move were the photometric inclinations used to derive T^i,/ values. Squares denote 
galaxies without bulges. Galaxies with evidence of mass discrepancies are marked by open symbols, and 
their Td / values are the maximum-disk values from Table 2. The horizontal errorbar in the lower right 
corner indicates the estimated error in the color adopted for each galaxy, (b) The same plot, but the different 
symbols denote whether a galaxy belongs to a clust^-|^or the field. Larger symbol sizes correspond to more 
luminous (and more metal rich) galaxies. The arrows show how points would move if their TuUy-Fisher 
distances were used instead of Hubble distances. 





Fig. 3. — Top left shows the projection of all Markov Chain Monte Carlo elements onto the po-Voo plane 
for the PI halo fit to ESO 322g82, to illustrate the strong degeneracy of these parameters. (Only 1 point in 
300 is plotted.) Top right shows the same projection of the 68% of points lying in the densest part of the 
parameter hypervolume. The other four panels show histograms of the same 68% of selected points projected 
onto the four parameter axes. The vertical line in each panel shows the best fit value. 
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Fig. 4. — Central density versus scale radius for the PI fits. The solid line is not a fit to the points, but 
shows the relation between the parameters for Voo = 200 km/s. Note a marginal tendency for the halos with 
higher central densities to fall below the line while lower density halos lie above it. 
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Fig. 5. — Concentration versus W200 for the NFW fits. Note how the uncertainties rapidly increase as the 
halos become less dense and more extended. The solid line illustrates the mean correlation found in the 
simulation of Bullock et al. (2001). 
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Fig. 6. — Slope parameter a versus vi for the power-law fits. The a values have an average value of 1.1, but 
have a wide range. 
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Fig. 7. — Fitted Td,i values vs. dereddened B — R when halos are included. The solid line in each panel 
is the best-fit correlation found by Bell et al. (2003). As in Figure 2(a), the solid squares denote galaxies 
without bulges. The error bars shown are the statistical errors presented in the tables, and do not include 
systematic errors. The horizontal error bar in the upper right corner of panel a) is the estimated color 
uncertainty for each galaxy. 
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Fig. 8. — Comparisons between the Td,i values fitted for the different halo models. Filled points mark 
galaxies which have '^d,i > 0.5 for both of the relevant models. Open circles indicate galaxies which do not 
require halos for either of the models. 
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Fig. 9. — Comparisons between the fractional halo masses for the different halo models. Open and filled 
points have the same meanings as in Figure 8. 
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Fig. 10. — Inner halo densities (A„/2) as a function of maximum halo velocity for the NFW halo function. 
Filled circles denote galaxies which have best-fit Td,! > 0.5, the open circle is for a galaxy for which no halo 
is required, and the points with error bars, but no circle are those for which T o.i < 0.5. The large error bars 
are for galaxies for which is poorly constrained by our data. The solid line is the predicted correlation 
from Zentner & Bullock (2002) and the error bar in the upper right hand corner is the A„/2 uncertainty 
estimated from Bullock et al. (2001). 
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Fig. 11. — Comparisons between velocity map rotation curves (solid error bars) and pseudo-slit rotation 
curves (dotted error bars) for 4 galaxies, (a) ESO 322g82 is an example of good agreement between the 
rotation curves, (b) ESO 322g36 has strong spirals which cause a difference between the approaching and 
receding sides of the pseudo-slit data, (c) ESO 438g08, the poor correspondence between the two curves 
arises because the photometric major axis is misaligned from the kinematic major axis, (d) ESO SOlgOl: 
the data along the major axis simply does not extend to large radius. 
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Fig. 13. — Comparisons of the dark matter halo masses for each of the halo models derived using pseudo- 
slit versus map data. Open circles indicate galaxies which do not require halos for either the pseudoslit or 
velocity map fit. 
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Table 1 

Palunas & Williams (2000) galaxies rejected for this study and reasons for exclusion. 



Galaxy Galaxy Galaxy 



Abell 1644d83^ 


ESO 


377gll3 


ESO 


444g21^ 


ESO 216g202 


ESO 


381g5li 


ESO 


444g86i 


ESQ 267g30i 


ESO 


382g58i'3 


ESO 


445g353 


ESO 269g6li 


ESO 


383g023 


ESO 


445g582'3 


ESO 322gl9i 


ESO 


383g883 


ESO 


445g8li 


ESO 322g48i 


ESO 


435g50i 


ESO 


SOlglli 


ESO 322g87i 


ESO 


436g39i 


ESO 


501gl52'3 


ESO 374g023 


ESO 


437g30i 


ESO 


501g862 


ESO 375gl23 


ESO 


437g342 


ESO 


509g9li 


ESO 375g29i 


ESO 


437g54i 


ESO 


510gll3 


ESO 376g02i 


ESO 


441g223 


ESO 


572gl73 


ESO 376gl0i'2'3 











^ inclination > 75° 
^ sparse Ha emission 
lack of central Ha emission 
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Table 2 

Best-fit Xbmin values, parameter values, and uncertainties for stars-only models. 



# 


Galaxy 


2 

X6,min 


Tr,,/(M0/L0) 




Ts,7(Mq/L0) 




1 

J. 




±.000\J 


9 "^9 


u.uy 


1 78 
±.10 


49 


9 




4:. 4:000 


9 1 S 


n 1 9 

u . xz 


9 18 


19 

U. 


Q 
o 






1 m 










T^.^n 9fi7o-9Q 


0.\J-L 1 U 


9 49 


1 RQ 
J. .uy 


9 49 


1 fiQ 
-L .uy 




p^n 91^80- '^7 


4 R79fi 


9 4'^ 


n fi8 






R 
U 




Zr.UOOU 


1 1^7 


n 07 


9 7'^ 


IQ 
u. ly 


7 


F.^n "^l 7o-41 


4 ^7RR 


9 47 


fl 04 


1 "^4 


0*^ 

U.UO 


s 

o 


F.SO 399p-'^fi/Nr;r! 4'i7'i 

J— /OVy KiLiLi^Kiyj j L \ V_i V_ ■ r±(J 1 tj 


4 Q4'^0 

'-t.ij'-t^yj 


1.40 


n 43 


1 in 


35 


q 




4 1 ^7fi 


1 '^Q 


n 0"^ 






1 fl 




9 49"^^ 


1 "^8 
1 .00 


nil 

u. 11 


1 "^8 
1 .00 


oil 

U. 11 


11 


ESQ S22ff45 


2.5188 


1.87 


0.01 


1.87 


0.01 




F,^n '?99cr7fi 


9 noi 9 


9 1 fi 


n 14 


1 1 


48 








0.00 


n 1 1 


1 8"^ 


7fi 
u. / u 


1 4 


F,Sn '^99o-89 /Nnr 4fi7Q 


4 07^0 


9 89 


n 1 n 


9 fn 

z . uo 


on 

u. yu 


i-O 




9 71 70 
Z. ( 1 ( u 


9 7'i 


n 77 


1 1 97 
± ± .z 1 


fi 48 


1 fi 
lu 


FCin '?9'^a-97 


9 S7'^9 


1 .yu 


19 
u. xz 


ii nfi 

o.uu 


07 
u.u / 


1 7 


F^n '^9'^ct'^Q 


1 9'^9 


1 "19 


1 00 








F^D '^9'^(r49 




1 117 
±.01 


09 


1 <\7 
±.0 1 


09 
u.uz 


1 Q 


F^D '^9'^o-7'^ 




1 "^Q 


97 






9fl 


F^D '^74o'n'^ 


491 n 


1 8n 


n OQ 

u.uy 






21 


Fsn 375p-09/Tr; 9559 


9 2998 


2.20 


20 


1 97 

± .Zi f 


58 


99 




O.U 1 


fi 1 7 


9 7^? 






9"? 




9 4409 


fi l^fi 
u . 


1 0*^ 
1 .uo 


4fi 


u.uu 


94 




1 1 ^^Q1 R 

± ± .Ot/ ±u 


9 


RR 
u.oo 


9 f^4 

Zi .0^ 


RR 
u.oo 






94Qfi 


9 Q9 
z .y z 


90 


9 Q9 

Zi .c/Zi 


90 


9R 




1 9074 


74 


1 47 


74 


1 47 


27 


ESQ 438g08 


2.0276 


0.92 


0.26 


2.35 


1.40 


28 


ESQ 438gl5 


2.7118 


1.67 


0.10 


2.02 


0.09 


29 


ESO 439gl8 


4.5172 


2.51 


0.35 


2.51 


0.35 


30 


ESQ 439g20 


2.0360 


3.33 


0.12 


3.33 


0.12 


31 


ESO 444g47 


2.0842 


3.27 


0.18 


3.27 


0.18 


32 


ESO 445gl5 


3.3342 


2.80 


0.19 


2.80 


0.19 


33 


ESO 445gl9/IC 4319 


2.4596 


2.28 


0.02 


2.28 


0.02 


34 


ESO 445g39/NGC 5298 


4.2358 


5.70 


0.68 


2.87 


0.78 


35 


ESO 446g01 


2.6448 


1.91 


0.46 






36 


ESO 501g01 


1.9834 


3.34 


0.72 


3.34 


0.72 


37 


ESO 501g68 


2.1836 


4.07 


0.27 


4.98 


0.17 


38 


ESO 502g02/NGC 3463 


1.9144 


2.37 


0.02 


2.37 


0.02 


39 


ESO 509g80/IC 4298 


3.6484 


3.73 


0.17 


2.92 


0.13 


40 


ESO 569gl7/NGC 3453 


1.6434 


2.26 


1.79 
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Table 2 — Continued 



# _ Galaxy x^.^in T,,,j(Mq/Lq) ax^ , TB_,7(Me/Le)_ _ arn., 
Maximum allowed T/ values from fits with the outer data discarded. 



10 


ESO 322g44 


2.1920 


0.98 


0.04 


0.98 


0.04 


19 


ESO 323g73 


1.1802 


1.31 


0.24 






25 


ESO 437g04 


3.2294 


2.47 


0.24 


2.47 


0.24 


31 


ESO 444g47 


1.2608 


2.67 


0.06 


2.67 


0.06 


35 


ESO 446g01 


1.1288 


1.67 


0.47 
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Table 3 

Best-fit 6^ values, parameter values, and uncertainties for PI models. 



Galaxy 





1D,/(.M0/L0j 




0j 


PO (M0/PC j 


Tc (kpc) 


Woo (km/s) 


1 


1.3822 


1.5e-06 ? ?'^ nR 


1.50 } 

-^•'-'^ 1.50 




21 o ?8 

0.12 


0.6 


176.46 


2 


4.0878 


0.95 

U.o4 


919 0.47 

0.00 




7.6e-02 ^-2^-02 

D.7e — 02 


4.0 }21-8 

1.5 


258.44 


3 


1.8988 


0.56 






n 1 1 

/.Ue — Uz 


2.0 

0. / 


152.39 


4 


5.3118 


1 60 S l^ 

j^.Kjyj 0.40 


1.95 ? If 

1.01 




3.1e-02 ?i'="Ho 

(j.j-v.. '-'■^ 1.8e— 02 


6.9 iV 
"."^ 3.9 


286.14 


5 


3.9596 


9.2e-07 S-?^ 






0.20 


0.1 


181.35 


6 


1.9840 


1.2e-08 ? 

i ,ze — Uo 


0.28 n OG 




9 09 2.48 

i . 04 


0.4 O;? 


151.97 


7 


4.5604* 


9.1e— 02 


1.43 

0.21 




1.9e-03 

1.2e— Uo 


1.2e+04 Jj^^ 


1.2e+05 


8 


4.4864 


0.23 23 


0.77 g:it 




n 1*^ 0.11 


2.2 o'g 


185.04 


9 


3.3828 


l.le-07 9-i9 ns 

o.7e— Oo 








4.1 ^1 

0.0 


161.18 


10 


1.5610 


0.42 S ?^ 


0.42 g-f^ 




1 91~ 

*J • 0^ n Qn 

U.OU 


0.9 


123.76 


11 


2.0248 


1.40 HI 

U. lo 


0.90 g-fg 




1.2e-02 ll'r?ii 

o.ze — UO 


2.6C+04 

z.oe-hu^i 


6.56+05 


12 


2.0442 


1.58 


1-69 °oM 

U.OZ 




6.0e-03 

o.oe — Uo 


15.2 


272.67 


13 


3.0918 


2.0e-05 


4.19 




0.27 g i 


2.5 fi 


305.55 


14 


4.0594 


9 oc 0.39 


9 77 0.79 
65 




0.20 lil 


0.9 li 


94.99 


15 


2.5246 


0.92 


0.92 




1.98 


rv 5.3e-02 
6.5e— 02 


207.34 


16 


2.6860 


8.0e-07 g:go_o7 


2.89 




n 1 7 0.53 

U. ± 1 Q le— 02 


2.5 


239.32 


17 


4.3368 


7.6e-08 9 i^ 08 

/ .D6 — UO 






4.0e-02'|-^^-"2 

±.oe — Uz 


2.9 g-f 

U.O 


136.76 


18 


1.5142 


0.78 8-i 


0.78 




8.OC-O2 ^-7-02 


1.8 


118.18 


19 


2.8006 


1 07 6.3e-02 

O.De — Uz 






2.0e-02 ?-2e-03 

1 .oe — Uo 


1.7e+04 


5.6e+05 


20 


3.1678 


1.6e-07 f ^ro7 






0.10 


1.8 g-? 


135.16 


21 


2.0276 


1.43 g }? 


2.62 9,ii 

0.47 




1.4C-02 ? n''"^? 

1.9e— Oo 


3.9e+04 ^^ttnl 

1 w i j.9e+u4 


l.le+06 


22 


2.2082 


5.9e-08 ff^_os 






w.tju g le — 02 


1 7 0-2 
1.1 0.2 


287.02 


23 


2.4352* 


3.79 ^ 


1.94 




0.11 i-ff °4 


1.5 ?'i 

1.0 


117.94 


24 


11.8492 


1.08 


1 51 1-41 




0.16 f8-02 

4:. oe — Uz 


1.9 "oi 

U.O 


174.89 


25 


2.4646 


1.7C-06 

1.7c — Od 


9.8C-07 Sj-3.|?_ 


-07 


n 90 3.6e — 02 
9.2e — 02 


^■'^ 0.2 


253.65 


26 


1.0336 


1.8c-08 ?'sl n« 
i .oe — Uo 


1.8e-08 yi^S 


-08 


0.19 3-?^-02 

4. ie — Uz 


1.5 01 


151.99 


27 


2.0276* 


Q QO 5.7e-02 
"•^"^ 6.6e-02 


9 qc 0.61 
.^•<J>' 0.63 




Q fip.n7 2 0e+09 


9.4e-02 ai_02 


0.01 


28 


2.5456 


"•^'J 0.32 


1 09 0.24 
-■-■^^ 0.18 




1 Kp no i.Oe-02 

i.Oe-UZ 8.6e-03 


9.8 


282.29 


29 


3.8966 


n 62 0-24 

U.OZ 0.21 


n 62 °-24 

U.OZ 0.21 




0.46 0:28 


1-7 g:l 


266.99 


30 


1.7666 


0.75 If, 


0.75 




0.80 


1 ri 0.2 
J-.U 0.1 


213.50 


31 


1.6208 


0.46 


0.46 Jii 




Q 8.4e-02 
'^•-'^-'^ 8.1e-02 


9 9 1.9 
^•^ 0.6 


168.11 


32 


3.0656 


1 70 0.29 
^■'^ 0.29 


1 70 0.29 
■■■•'o 0.29 




14 O '*^ 

^•-■■^ 0.11 


1 c 1.8 
0.8 


144.76 


33 


2.2896 


2.9e-06 Ul-06 


1.00 0.72 




Q iQ 5.0e-02 
^•■•■^ 5.1e-02 


9 4 0.4 
0.3 


241.52 


34 


4.1728 


3.3C-06 i:ito6 


3.29 lit 




"•^^ 0.12 


3.7 \i 


404.20 


35 


2.2198 


. 44 5.9e-Q2 
6.0e-02 






l.le-02 3 3e_03 


9.7 il 


233.82 


36 


2.2832 


7.2e-10 2!9e-18 


6.7e-07 


-07 


O.ae UZ 4 7e_02 


9 r, 0.5 
0.2 


136.94 


37 


2.1342 


2.40 \il 


2.79 2-19 




4.99 


0.2 J«2' 


117.32 


38 


1.6984 


0.63 g:i 


0.63 °il 




9 10 1.99 

1.63 


0.6 8:f 


191.69 


39 


3.4298 


0.41 'o.U 


0.41 




97 4.9e-02 
5.3e-02 


9 9 0.2 
^.^ 0.1 


265.12 


40 


1.4480 


6.1e-10 0:?o_io 






8.24 HI 


f^ 4 4.9e-02 
5.4e-02 


243.78 



36 



^Vca is not a fitted parameter. It has been calculated using the best-fit po and rc values to facilitate 
comparison with the other halos' velocity parameters. 

*Indicates that the halo fit is no better than the stars-only fit at the 95% level. 
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Table 4 

Best-fit 6^ values, parameter values, and uncertainties for NFW models. 



Galaxy 






Tr rCMn/Ln) 


c 


7.?onn 1 km /s ) 




1 


1.3946 


3.9e-07 0:1^07 


1 13 


15.6 Itl 


107 Q 29.7 


10.9 


2 


4.1132 


'-■■^'J 37 


1 0.39 

l.OO g 4.3 


6.0e-02 ?-23_p2 


3.8e-t-03 3 6e+o3 


8.5e+04 


3 


1.9000 


n 27 




17.9 Ig5| 


120.0 o 


8.9 


4 


5.3142 


1-10 gii 


1-97 8:?l 


0.15 


4.0e+03 ^ o^+o^ 


3.5e+04 


5 


4.0182 


5.0e-08 lll_os 




16 7 2-43 
J-".' 3.06 


122.0 fl 


9.8 


6 


1.9934 


0.29 


0.29 


44.3 i«46'' 


69.6 ^4^7 


2.1 


7 


4.6242* 


1 01 0.39 

l.Oi g 24 


1 oi 0.39 
l.Ol 24 


ill 


151.9 .^i 


17.7 


8 


4.4812 


U.iZ g 


19 "•4'' 

'-'•-'-^ 0.12 


18.0 7;7g 


121.9 fd 


9.0 


9 


3.5878 


1 no 06 0-20 

i.UC-UO i.Oe-06 




4.06 


210.9 i^iyO 


69.3 


10 


1.5704 


n OO 5.2e-02 


n 80 5.2e-02 
U.SZ Q 26 


4.05 12^02 


145.5 


47.9 


11 


2.0276 


1.4e-06 


1.4e-05 ijo.os 


12.9 ifl 


172.2 ly^f 03 


17.8 


12 


2.0502 


^•'^^ 29 


1.58 g'ly 


5.7C-04 2-52_g^ 


1.6C+03 r2°,%3 


3.7C+06 


13 


3.0854 


1.6e-06 \£_oe 


2.86 }:°? 


10.8 


qro 3.0e+03 
OdO.O 223 1 


44.4 


14 


4.0478 


1-45 


1.42 l 


30.3 l?-^i 


91.5 \n 


4.0 


15 


2.4916 


0.82 " i 


0.82 0-26 


39.8 til 


112.7 


3.8 


16 


2.6952 


4.5e-07 O'ltgy 


1-71 Ifs 


17 Q 6.35 
-'- ' 5 46 


160.1 ill 


12.3 


17 


4.3938 


l.Oe-07 "i::^^ 




5.87 1°^ 


140.2 i26g4 


31.8 


18 


1.5374 


0.97 0-21 


0.97 gji 


8.75 til 


84.9 2|2 


12.9 


19 


2.8008 


0.62 1-%-"^ 




0.17 g fi 


3.5gH~03 QD_Ln^ 


2.7e+04 


20 


3.2332 


4.8e-07 ^i^^I^^ 




11.6 H? 


102.5 ii 


11.8 


21 


2.0202 


1 1 0-34 

^ 4 5e — 02 


n 1 1 0-34 

^ 4.5e— 02 


18 q 2-17 

6 60 


112.527-9 


7.9 


22 


2.1980 


1.2e-06?fo_g6 




1 3!80 


211.4 III 


15.5 


23 


2.4402* 


6.56 


3.46 


2.06 iil 


2.6e-03 182,0 g3 


0.0 


24 


11.8384 


0.58 8 i 


0.58 0-28 


19.8 2 I0 


127.2 


8.6 


25 


2.5360 


1.5e-06 ?i3_gg 


1.5e-06 o-itoe 


12.8 lil 


216.6 29-2' 


22.5 


26 


1.0658 


^•■■■"^ 1.29 


1 c 0.15 
Z.IO ^ 29 


Qr 10.31 


267.6 


1.56+02 


27 


2.0276* 


Q9 6.3e — 02 
^■^^ 6.4e-02 


0.41 
0.70 


171.4 0i2i 


8.0C-02 8 5e_g2 


0.0 


28 


2.5472 


7.2e-02 7;2e-o2 


1 29 " 
J--^^ 0.20 


1.8e-04 fii_04 


3.1e+03 2.8e+03 


2.3e+07 


29 


3.9182 


2.8e-06 if,_oe 


2.8e-06 ^rLoe 


QR Q 4.18 
<J"."5 6.87 


157.1 ti 


5.8 


30 


1.7894 


5.1C-07 0:63_o7 




35 2.26 
•^^^-S 5.95 


127.8 li 


4.8 


31 


1.6606 


1 ■^fi 0-31 




4.01 i:?f 


248.6 350-8 


82.7 


32 


3.0718 


1 7<^ 0-21 
-■■•'3 0.34 


1 7'^ 0-21 
J-.' 3 0.34 


1/1 19.36 
9.90 


101.8 94-2 


9.2 


33 


2.2710 


1 -26-07 i;2e-07 


9.8C-02 9:ie_02 


Ifi 9 2-62 
J^O.^ 4.00 


167.6 27.4 


13.8 


34 


4.2076* 


1.4C-04 


o 04^ 0.66 
^•^^ 0.75 


8.8C-06 8:8eIo6 


1 ^p-l-04 3.8e+03 
i.0e+U4 1 2e+04 


2.3C+09 


35 


2.2232 


T .,1 4.5e-02 

i.Oi g ]^2 




1.86 f i 


47f; 775.0 
*'0.Z 318.4 


3.4C+02 


36 


2.3106 


fi 9r 07 0-43 
D.ZC U/ 2.0e-07 


6.2C-07 ^:43_„, 


11 3 i -''5 

^^•'J 1.80 


103.1 fl 


12.2 


37 


2.1400 


9 4-? 1-91 
0.98 


4.28 


t- 07 4.8e+02 


1 C4 152.0 


38.3 


38 


1.7024 


0.37 


0.37 


AO -i 10.78 


106.7 li 


3.4 


39 


3.4314 


2.0C-06 ;:30_o6 


2.0e-06 ;;30_og 


20.1 lil 


179.2 fl 


11.9 


40 


1.4178 


9.0e-09 0:30_g9 




69.4 


105.5 f ^ 


2.0 



38 



^Rg is not a fitted parameter. It has been calculated using the best-fit V200 and c values to facilitate 
comparison with the PI scale length parameter. 

*Indicates that the halo fit is no better than the stars-only fit at the 95% level. 
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Table 5 

Best-fit 6^ values, parameter values, and uncertainties for POWER models. 



Galaxy 





1 d,i(Mq/Lq) 




vo (km/s) 


a 


1 


1.4170 


9.3e-08 ^i'i 

i.oe— ib 


U.oo 


158.2 ^if. 


1 2 °i 

0.2 


2 


4.1074 


1.47 "-50 


1.79 o « 

U.OO 


162.5 §o g 

00.0 


0.7 H 

0.7 


3 


1.9188 


0.56 ^1? 




134.2 19-9 

" 25.5 


0.4 


4 


5.3118 


1 .56 Q 


1.56 n 


161.8 li'i 


0.8 


5 


4.0946 


0.89 0^22 




151.4 ?i57 


1 0.1 

2 


6 


2.0340* 


1-41 °Al 


U.4D 


48.7 


1.6 

1 .b 


7 


4.6356* 


1 .T5 Q 22 


1 .T5 Q 22 


140.7 Pi 


1 1 °-2 

U.O 


8 


4.4868 


0.43 0^11 


0.43 


156.5 lit 


1 /I 7.1e-02 


9 


3.5878 


0.39 0-21 




117.1 1?:? 


0.8 Sj'i 
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Table 5 — Continued 



Galaxy 6^ TdAMq/Lq) Tb,i{^qM (km/s) a 



*Indicates that the halo fit is no better than the stars-only fit at the 
95% level. 



